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NUMERICAL  SIMULATIONS  OF  MAGNETO-HYDRODYNAMIC  CHEMICALLY 
REACTIVE  FLOWS  IN  THREE-DIMENSIONS  ON  A  VECTOR  COMPUTER 


I.  Introduction 

Large  numerical  simulations  of  three  dimensional  flows  have  been 
performed  at  the  Naval  Research  Laboratory.  The  limitations  of  the  present 
day  computers  with  a  word  storage  of  about  10^  words  and  a  speed  of  about  10^ 
operations/sec  impose  severe  restrictions  on  the  kind  of  problems  one  can  do 
in  a  reasonable  time  on  a  computer. 

The  storage  restriction  can  be  somewhat  alleviated  by  using  discs  as  an 
Intermediate  storage.  But  still  relative  coarse  resolution  of  15-30 
meshpoints  per  dimension  is  a  serious  drawback.  It  can  be  partially  overcome 
by  using  an  appropriate  coordinate  system  and  a  closely  connected  difference 
scheme. 

The  requirements  on  a  good  numerical  scheme  in  multidimension  can  be 
stated  as  follows: 

1)  at  least  second  order  in  space  and  preferably  also  in  time 

2)  as  non-  diffusive  as  possible 

3)  the  convection  terms  have  to  be  monotonic.  This  implies  the  use  of  a 
(nonlinear)  hybridization  with  a  first  order  scheme 

4)  adjusted  coordinate  system,  stretched  mesh 

5)  all  points  should  be  closely  connected.  This  excludes  in  3-D  the  use 
of  a  leapfrog  scheme. 

6)  the  scheme  should  be  non-time  split 

7)  the  scheme  should  be  conservative  to  the  extent  possible. 

Many  problems  require  a  transport  of  many  chemical  species  which  react 
chemically  with  each  other  at  one  raeshpoint.  The  resulting  ordinary 
differential  equations  are  generally  stiff.  The  scheme  used  here  is  a  second 
order  predictor,  corrector  method. 


Manuscript  aubmitted  February  8,  1982. 


1 


In  chapter  II  the  hydrodynamic  scheme  will  be  briefly  explained  with  a  1- 
D  example.  Chapter  III  shows  the  implementation  of  the  scheme  in  a  3-D 
code.  Specially  the  problem  of  keeping  the  magnetic  fields  divergence  free 
and  the  problem  of  high  characteristic  speeds  are  discussed.  Chapter  IV 
presents  the  numerical  scheme  for  the  ordinary  differential  equations  and 
shows  the  data  flow.  Chapter  V  gives  an  overall  flowchart  for  a  multifluid 
code  and  for  interweaving  computations,  writing  and  reading  discs.  Chapter  VI 
shows  the  results  of  3-D  computations  for  the  plasma  flow  around  the  earth. 

The  example  demonstrates  clearly  that  even  with  a  coarse  resolution  it  is 
possible  to  get  detailed  results. 


II.  One  dimensional  example 

The  numerical  scheme  (PHOENIX)  presented  here  Is  the  simplest  possible 
second  order  scheme.  The  hydrodynamic  variables,  density  and  pressure  are 
defined  at  oeshpoints  and  time  t,  the  velocities  as  defined  at  midpoints  and 
time  t  -  ,  the  magnetic  fields  are  given  at  the  same  points  In  space  as 

the  velocities,  but  at  time  t.  PHOENIX  Is  an  explicit  code,  therefore  the 
time  step  limitation  is  given  by  the  Courant  condition. 

The  scheme  will  be  explained  on  the  one  dimensional  continuity  equation  for 
the  number  density  In  cartesian  coordinates 


3n 

3t 


Let  us  define  the  following  operators 


"  ■  i  (f.i +  f.i  - 1> 


,  ".i*  ‘  Vi  •  fj 


5  <V 


pdm  (v,  f)  »  (m  -  n  sign  (v)  7  )f 


Tne  numerical  scheme  for  solving  the 

n  “  n  “  %x  *(v  Pdm  Cv.n)) 
where  ,,AM  denotes  the  value  at  t  +  dt. 

The  POM  operator,  the  so  called  partial  donor  cell  operator,  will  now  be 
-*xplatned.  The  diffusion  coefficient  n  has  to  he  chosen  such  that  n  remains 
monotonic.  This  means  that  the  convected  n^  has  to  be  bounded  by  n^+j,  and 
nj_l*  006  can  easily  see  that  n  -  0  gives  no  diffusion  and  n  ■  1/2  the  donor 
cell  method.  The  partial  donor  cell  method  [1]  determines  n  nonllnearly.  One 
should  remark  here  that  this  method  gives  no  clipping  for  v  -  0  ,  In  contrasts 
to  flux  corrected  transport  [2,3]  FCT.  The  method  can  be  easily  used  to 
hybridise  higher  order  schemes  [4],  In  the  PHOENIX  scheme  the  formula 
for  n  7f  is 


fi  '  f.i  i  i 

where  n  is  the  diffusion  coefficient, 
continuity  equation  can  be  written  as 

A  At  4- 
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f|Vf  «■  —  sign  (Vf)  max  |0, |Vf| 

-  6/4^sign  (Vf+)+  sign  (Vf)^(l  +  sign(v))  |Vf+| 

-  6/4  ^sign  (Vf-)  +  sign  (Vf)^  (1  -  sign  (v))  |Vf-j J 

where  6  =  .6  is  usually  used  to  maintain  a  small  amount  of  diffusions. 

(6=o  the  donor  cell).  One  may  also  note  that  an  extremum,  the  maximum 
diffusion,  is  taken  in  order  to  avoid  overshoots. 

PHOFINIX  uses  velocities  and  pressure  as  variables.  The  momentum  and  energy 
equation  are  solved  in  the  same  fashion  as  explained  above.  It  is  easy  to 
conserve  total  mass  numerically.  All  other  physical  quantities  like 
momentum  or  total  energy  are  conserved  to  the  order  of  the  scheme.  Energy 
conservation  is  used  as  a  check.  Theoretically  the  nonconservation  is  small, 
hut  not  necessarily  in  practice  because  of  the  occurrence  of  discontinuities. 

In  general,  it  is  not  advisable  to  conserve  total  energy 
numerically  because  the  error  in  computing  the  pressure  (temperature)  can 
become  very  large.  The  temperature  is  a  critical  variable  for  reactions  and 
radiation  and  has  to  be  as  accurate  as  possible.  A  scheme  which  does  not 
conserve  total  energy  numerically  has  to  loose  energy  in  order  to  be  stable. 
In  order  to  show  that  the  code  works  properly  in  one  dimension,  the  numerical 
solution  of  a  diaphram  problem  as  posed  by  Sod  [51  has  been  used  as  a  test. 
The  parameters  are  as  follows 

Pj  =1.  P2  =  .1 


Y 1  =  1.4 


P2  -  .125 

V2  *  0 

Y2  -  1.4 


The  resulting  machnumber  is  1.62.  Figures  1,2,3  show  the  densities,  presures 
and  velocities  after  the  shock  traveled  a  quarter  of  the  mesh  (t=.147  sec). 

One  solid  line  gives  the  analytic  solution  of  the  Riemann  problem,  the  other 
solid  line  shows  the  results  of  a  high  resolution  calculation  with  180  meshpoints, 
the  dashed  the  results  for  25  meshpoints.  One  can  see  that  the  results  for  the 
high  resolution  case  are  quite  good.  For  the  low  resolution  case  the  main 
features  are  still  represented.  This  is  the  kind  of  resolution  one  will  find 
in  a  typical  3-D  calculation. 


T  — 
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Diaphragma  at  .147  sec:  pressure  profiles 
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III.  Extension  to  two  and  three  dimensions 

The  PHOENIX  scheme  can  be  easily  extended  to  two  and  three  dimensions, 
but  several  points  should  be  made. 

Shockwidths  should  be  constant  in  the  number  of  points  over  which  the 
shockfront  is  computed.  The  shock  should  not  be  more  than  4  meshpoints 
wide.  In  general  the  meshspacing  is  not  the  same  for  different  directions, 
therefore,  it  follows  that,  the  artificial  shock  pressure  used  in  the  code  is 
not  isotropic  in  physical  space.  PHOENIX  gives  a  shockwidth  of  2-3 
meshspacing.  In  a  magnetohydrodynamic  flow  the  magnetic  field  has  to  remain 
numerically  divergence  free.  This  can  be  achieved  in  two  ways,  either  one 
computes  the  vector  potential  or  one  uses  the  curl  of  the  electric  fields  to 
advance  B.  In  many  2-D  problems  the  first  approach  is  useful  because  there  Is 
only  one  component  of  the  vector  potential.  In  3-0  the  latter  approach  is 
more  suitable.  The  PHOENIX  scheme  computes  the  electric  fields  in  the 
magnetic  field  calculation.  This  approach  avoids  the  taking  of  second 
derivatives  which  may  introduce  numerical  noise.  In  the  compuatton  of  the 
electric  field  the  PDM  operator  is  used.  It  results  in  a  resisitivity  which 
causes  B  to  be  monotonic  (compare  [2]). 

In  most  real  MHD  problems  there  exists  regions  where  the  Alfven  speed 
2 

ca  =  8  /^P  is  high,  therefore  requiring  very  small  timesteps.  The  time  step 

can  be  increased  In  two  ways,  either  by  using  an  implicit  scheme  for  computing 

B  or  by  limiting  the  characteristic  velocity  to  a  velocity  c,.  The  first 

",1m 

approach  is  feasible  in  2-D.  It  Involves  solving  impltcity  a  diffusion 

2 

equation  with  an  diffusion  coefficient  c  dt  .  The  second  approach  consists 

d  2  2 

of  increasing  the  density  perpendicular  to  B  by  a  factor  1  + 

thus  decreasing  the  accelerations  perpendicular  to  the  magnetic  field.  The 

limiting  velocity  c? ^  should  be  larger  than  any  physical  velocity  involved. 
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IV.  Chemical  reactions 

The  ordinary  differential  equations  resulting  from  chemical  reactions  at 
each  meshpoint  can  be  written  as 


dy 

dt 


-  ay  4-  b 


where  y  represents  the  species  array.  The  coefficient  a  and  b  depends  also  on 
time  via  species  concentrations.  Assuming  a,  b  to  be  constant  during  a 
times tep  t,  the  exact  solution  is 


y  -  [y„  +  h  (.“  -  Dl 

The  numerical  approximation  to  this  exact  solution  was  chosen  such  that  it 
will  give  the  correct  asymptotic  solution 


y  *  -g-  for  at »  1 

In  order  to  acchieve  this  eat  is  approximated  by 

at  1  2  2 

e  »  1  +  at  +  ■j  a  t 

(Some  integration  schemes  use 

eat  »  (1  +  -|at)/(l  -  ^  at) 
which  will  lead  to  ocillations  for  at  »  l) 

As  a,  b  are  time  dependent  one  has  to  use  a  predictor  and  corrector 
step. 


Predictor  step 

-  <y0  +  bo  °  + 1  v)t)/(l  +  V  +  i  %  t2) 

Corrector  step 

yi  "  (yo  +  1  (bo  +  be  (1  +  V))t)/(1  +  I  (ao  +  ae(1  +  aet^)t) 
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The  difference  betwen  ye  and  y*  determines  the  time  step. 

The  setting  up  of  the  differential  equations  and  the  coding  for  the 
numerical  Integration  Is  achieved  by  an  interpretive  program  ARIS  [6J.  The 
Input  to  this  program  are  the  chemical  reactions  and  reaction  energies.  The 
resulting  FORTRAN  code  is  completely  vectorized.  Decisions  are  made  by  the 
use  of  sign  functions  rather  than  IF  statements.  This  permits  great 
improvements  in  speed  as  the  sign  iunctton  statements  are  vectorized. 

The  chemical  reactions  have  to  be  solved  at  every  meshpoint  for  a  certain 
time  interval  tc^e  •  As  the  species  concentrations  and  temperature  are 
different  at  every  meshpoint  the  number  of  timesteps  needed  for  a  given  error 
varies  from  point  to  point.  As  the  computations  are  fairly  complicated,  the 
number  of  auxiliary  arrays  Is  large,  therefore  the  simultaneous  solutions  are 
set  up  for  an  array  of  length  m  <  n  where  n  is  the  number  of  meshpoints. 
Usually  m  =■  200.  The  data  flow  is  shown  in  the  following  flowchart 
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In  this  way  the  chem  array  is  always  filled  until  the  last  point  is 
completed.  This  optimization  achieves  an  efficiency  of  about  75 t  of  the 
maximum  speed  (about  3.10^  op/ sec  at  the  ASC  of  Texas  Instruments). 

V.  Flow  charts  and  data  flow 

The  following  flowchart  gives  the  general  computational  scheme  for 
solving  the  magnetohydronamic  equations  for  one  timestep.  As  the  PRORNIX 
scheme  is  explicit,  the  timestep  has  to  be  determined  by  the  Courant 
condition.  Furthermore,  the  stability  conditions  require  that  first  all 
p,  B  computed,  afterwards  the  velocities  are  updated.  In  addition,  a 
limiting  velocity  c  is  used  here 
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In  a  multifluid  code  each  of  the  fluids  can  have  their  own  time step.  They 
interact  with  each  other  through  collision  terms.  The  chemistry  is  called 
Independently. 

Some  of  our  large  simulations  require  a  data  base  which  exceeds  the  available 
core  storage  of  about  10^  words.  This  requires  the  use  of  disks  as  an 
intermediate  storage.  The  computation  is  done  for  one  plane  in  the  third 
dimension  at  one  time. 

PDM  (or  FCT)  requires  the  code  to  keep  at  least  5  planes  in  core. 
PHOENIX,  because  of  the  definition  of  v^  at  midpoints,  needs  one  more  plane. 
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The  reading  and  writing  has  to  be  done  simultaneously  with  the  computations. 
This  adds  two  more  planes.  Therefore  ft  planes  have  to  he  in  core.  The 
restriction  on  the  number  of  points  nj  x  nj  is  therefore  given  by 

8  x  T.  x  nixn-><  core  -  auxilary  variables 
var  1  4 

where  L  Is  the  total  number  of  variables  used, 
var 

PDM  (or  FCT)  require  the  computation  of  higher  order  differences.  For 
the  third  direction  these  computations  have  to  he  repeated.  This  Increases 
the  number  of  computations  by  about  102. 

There  are  about  ISO  operations  per  point  and  variable.  The  time  per 

-8 

operation  on  the  ASO  of  Texas  Instruments  Is  about  3  x  10  sec, 

,6 

therefore  the  computing  time  per  point  and  variable  Is  “  5  10  sec.  Assume 

n,  x  n-  *  650,  and  the  number  of  variables  to  be  30,  such  that 
\  i  4 

j l  x  ni  x  n?  »  2  x  10  .  This  gives  "  .1  sec  per  plane.  For  SO  planes  this 

var  i  c  . 

will  give  5  sec  per  tlmestep. 

The  1/0  operations  are  done  In  one  1/0  instruction  without  using 
additional  buffer  storage.  The  correctness  of  the  data  transfer  is  checked 
after  completion. 

The  core  storage  Is  used  In  a  revolving  fashion.  If  L  Is  the  Index  to 
reference  a  2-Ti  plane,  then  the  2-d  planes  of  all  variables  for  -L-  are  stored 
as  one  sequential  block.  The  start  Index  for  such  a  block  is  given  by 

T,f(h)  -  T.  x  nix  n?  x  mod(L  +  7,81  +  1 
va  r  1  '■ 

The  address  of  one  special  variable  n.  Is  then  given  by 

L  (I.)  ■*  Lf (L)  +  (n  -  1  )  x  ni  x  n2 
varn  1  c 

To  be  able  to  write  a  comprehensible  program  the  names  of  physical 
variables  are  set  with  equivalences. 

In  order  to  save  1/0  operations,  the  hydrodynamic  calculation  can  run 
forward  or  backward  in  the  third  dimension.  These  are  always  the  first  or 
last  three  planes  in  core.  The  following  flow  chart  will  show  the  data  flow 
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VI.  3-D  example 

As  an  illustration  of  the  capabilities  of  the  3-D  MHD  code,  the 
interaction  of  the  solarwind  with  the  earth  dipole  field  has  been  computed.  t 
[8]  The  earth  dipole  field  is  formed  by  placing  a  current  loop  in  one  of  the 
cells.  A  solar  wind  is  then  introduced  at  the  left  x  boundary.  It  has  a 
velocity  of  400  km/sec  and  carries  a  density  of  5  particles/cm' .  In  the 
example  shown  it  carries  also  a  southward  field  of  -2y  .  The  length  unit  is 
one  earthradius.  The  computation  shown  required  about  half  an  hour  of 
computing  time  and  simulated  two  hours  in  real  time  (about  2000  time  steps). 

The  mesh  is  2D  x  21  x  21  “  13000  points.  Several  difficulties  arose  durtng 
the  runs.  The  main  one  involved  the  internal  boundary  at  the  earth.  This 
problem  produced  large  currents  due  to  the  sharp  field  curvature  and  low 
resolution.  This  was  solved  by  massloading  the  cells  around  the  earth,  and  not 
allowing  the  mass  to  move. 

One  other  interesting  phenomena  occurred.  The  plasma  sheet  behind  the 
earth  contains  opposite  fields.  The  configuration  is  physically  unstable  to 
resistive  instabilities,  resulting  in  reconnection  of  field  lines.  Although 
the  code  Is  nearly  diffusion  free,  this  occurred  during  test  runs  showing  up 
as  reconnections  of  the  magnetic  field  lines,  a  result  that  is  thought  to 
cause  magnetic  substorms  when  physically  initiated,  [8] 

The  Figs.  4  and  5  show  density  and  pressure  in  the  noon-midnight  meridional 
plane.  The  tines  display  the  actual  grid  used.  One  can  see  the  density 
increase  at  the  bowshock  by  abdut  a  factor  of  3.  The  raachnumber  is  about  3. 

Also  one  can  see  a  slight  Increase  in  density  in  the  plasma  sheet  which  is 
about  about  one  cell  wide.  The  pressure  graph  shows  the  structure  a  bit 
clearer.  It  also  shows  the  side  lobes  of  higher  pressure. 

These  results  show  clearly  that  one  can  keep  detailed  structures  even 
with  a  relative  coarse  resolution.  They  also  demonstrate  the  almost  complete 
lack  of  diffusion. 
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Fig,  5  —  Bow  shock  with  IMF  *  -2y  ■  pr< 


VII.  Conclusion 

It  is  possible  to  perforin  realistic  3-D  simulations  with  present  day 
computers. 

The  use  of  an  adjusted  coordinate  system  in  combination  with  a  numerical 
scheme  which  closely  connects  neighboring  points  can  give  detailed  structures 
even  with  a  relative  coarse  resolution.  The  non  diffusivity  will  keep  those 
structures  during  a  calculation.  The  use  of  virtual  disc  storage  can  al  least 
partially  overcome  the  limitations  of  the  core  storage.  The  optimization 
achieved  by  the  vectorization  of  the  chemical  reactions  makes  it  possible  to 
do  realistic  3-D  flow  problems  which  involve  many  species  and  their 
interactions. 
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01CY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORRE YANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

01CY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAURAL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEMRA 
01CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  00) 

SANTA  BARBARA,  CA  93102 
01CY  ATTN  DAS I AC 
01CY  ATTN  TIM  STEPHANS 
0  ICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
01CY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIF.GO ,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSI1.ES  6  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  <>4088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 


LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01CY  ATTN  MARTIN  WALT  DEPT  52-12 
01CY  ATTN  W.L.  IMHOF  DEPT  52-12 
01CY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
01CY  ATTN  J.B.  CLADIS  DEPT  52-12 

LOCKHEED  MISSILE  &  SPACE  CO.,  INC. 
HUNTSVILLE  RESEARCH  &  ENGR.  CTR. 

4800  BRADFORD  DRIVE 
HUNTSVILLE,  AL  35807 
ATTN  DALE  H.  DIVIS 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
01CY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  TABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

01CY  ATTN  DAVID  M.  TOWLE 
01CY  ATTN  P.  WALDRON 
01CY  ATTN  L.  LOUOHLIN 
01CY  ATTN  D.  CLARK 

MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CY  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOIJLE 
01CY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
01CY  ATTN  R.W.  HALPRIN 
01CY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 


MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
01CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01CY  ATTN  W.  HALL 
01CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD. 

SANTA  MONICA,  CA  90404 
01CY  ATTN  E.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
31R  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

01CY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 


01CY 

ATTN 

P.  FISCHER 

R  &  D  ASSOCIATES 

01CY 

ATTN 

W.F.  CREVIER 

P.O.  BOX  9695 

01CY 

ATTN 

STEVEN  L.  GUTSCHE 

MARINA 

DEL  REY,  CA  90291 

01CY 

ATTN 

D.  SAPPFNFIELD 

01CY 

ATTN  FORREST  GILMORE 

01CY 

ATTN 

R.  BOGUSCH 

01CY 

ATTN  BRYAN  GABBARD 

01CY 

ATTN 

R.  HENDRICK 

01CY 

ATTN  WILLIAM  B.  WRIGHT,  JR. 

01CY 

ATTN 

RALPH  KILB 

OlCY 

ATTN  ROBERT  F.  LELEVIER 

01CY 

ATTN 

DAVE  SOWLE 

01CY 

ATTN  WILLIAM  J.  KARZAS 

01CY 

ATTN 

F.  FAJEN 

OlCY 

ATTN  H.  ORY 

01CY 

ATTN 

M.  SCHEIBE 

OlCY 

ATTN  C.  MACDONALD 

01CY 

ATTN 

CONRAD  L.  LONGMIRE 

OlCY 

ATTN  R.  TURCO 

01CY 

ATTN 

WARREN  A.  SCHLUETER 
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RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  FND  AVENUE 
NEW  YORK,  NY  10023 

OICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.O.  BOX  2351 
LA  JOLLA,  CA  82038 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  E.  FRIEMAN 
OICY  ATTN  E.A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  COODRIDCF  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

OICY  ATTN  DONALD  NEILSON 
OICY  ATTN  ALAN  BURNS 
OICY  ATTN  0.  SMITH 
OICY  ATTN  L.L.  COBB 
OICY  ATTN  DAVID  A.  JOHNSON 
OICY  ATTN  WALTER  G.  CHESNUT 
OICY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  M.  BARON 
OICY  ATTN  RAY  L.  LEADABRAND 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  J.  PETERSON 
OICY  ATTN  R.  HAKE,  JR. 

OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 

STEWART  RADIANCE  LABORATORY 
UTAH  STATE  UNIVERSITY 
1  DE  ANGELO  DRIVE 
BEDFORD,  MA  01730 
OICY  ATTN  J.  HLWICK 


TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 
OICY  ATTN  W.P.  BOQUIST 

TRW  DEFENSE  A  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDT 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 
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